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Abstract.  In  this  paper,  we  clarify  several  outstanding  issues 
concerning  the  predominant  mechanism  of  vorticity  gener¬ 
ation  in  mesoscale  convective  vortices  (MCVs)  in  weak  to 
modest  baroclinic  environments  with  nonzero  Coriolis  pa¬ 
rameter.  We  examine  also  the  corresponding  diabatic  heat¬ 
ing  profiles  of  the  convective  and  stratiform  components  of 
the  MCS  and  their  effects  on  the  concentration  and  dilution 
ofPV  substance. 


1  Introduction 

In  their  pioneering  paper  outlining  a  theory  for  the  mainte¬ 
nance  of  long-lived  mesoscale  convective  systems  (MCS)s, 
Raymond  and  Jiang  (1990)  stated  that  “the  convective  and 
anvil  regions  have  the  same  qualitative  effect  on  potential 
vorticity  distributions,  i.e.  potential  vorticity  in  the  lower  half 
of  the  troposphere  is  increased,  whereas  potential  vorticity 
near  the  tropopause  is  decreased”.  With  regard  to  the  ab¬ 
solute  vorticity  generation  in  mesoscale  convective  vortices 
(MCVs),  it  is  commonly  thought  that,  because  the  vortex  re¬ 
sides  in  an  area  of  stratiform  precipitation,  the  stratiform  pro¬ 
cesses  must  be  largely  responsible  for  generating  its  absolute 
vorticity  (Hertenstein  and  Schubert,  1991;  Johnson  and  Bar¬ 
tels,  1992;  Fritsch  et  ah,  1994;  Davis  and  Trier,  2002). 

“This  large  derivative  of  the  heating,  coupled  with  the 
longer  influence  time  associated  with  the  width  of  the  strat¬ 
iform  region,  allows  the  potential  vorticity  signature  of  the 
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stratiform  region  to  dominate  over  the  signature  of  the  con¬ 
vective  line”  (Hertenstein  and  Schubert,  1991). 

“Intensification  of  the  MCV  began  overnight  when  a 
lower-tropospheric  mesoscale  vortex  formed  on  the  north¬ 
ern  end  of  the  north-south-oriented  convective  line.  Inten¬ 
sification  of  the  midtropospheric  vortex  followed,  occurring 
in  response  to  the  development  of  a  stratiform  precipitation 
region”  (Davis  and  Trier,  2002). 

A  number  of  these  studies  have  indeed  provided  partial  an¬ 
swers  to  the  questions  regarding  the  generation  mechanism 
of  absolute  vorticity  in  the  MCV.  However,  as  Raymond  and 
Jiang  stated,  “further  work  is  needed  in  determining  the  con¬ 
tribution  from  each  effect  in  a  variety  of  mesoscale  systems.” 
The  purpose  of  this  paper  is  to  take  a  more  detailed  look  at 
the  vorticity  generation  mechanisms,  both  in  terms  of  PV  and 
absolute  vorticity.  In  particular,  we  will  take  a  look  at  the 
relative  contributions  to  the  PV  concentration  from  the  con¬ 
vective  and  stratiform  regions  of  two  simulated  MCSs. 

In  terms  of  the  absolute  vorticity  budget,  the  primary  in¬ 
tensification  mechanism  is  by  means  of  the  thermally  direct 
circulation  driven  by  diabatic  heating  within  the  convective 
system  and  the  associated  flux  convergence  of  planetary  and 
relative  vorticity.  When  referring  more  directly  to  potential 
vorticity  (PV),  the  intensification  can  be  explained  (at  least 
partially,  when  viewed  in  geometric  coordinates,  using  the 
material  form  of  the  PV  equation)  in  terms  of  the  vertical 
gradient  of  diabatic  heating  within  the  convective  system.  In 
either  case,  diabatic  heating  is  required,  whether  it  be  by  driv¬ 
ing  the  thermally  direct  circulation  or  producing  PV  by  ver¬ 
tical  gradients. 

The  aim  of  this  paper  is  to  resolve  a  number  of  ambiguous 
conclusions  regarding  the  sources  of  PV  in  MCVs.  In  their 
satellite-based  study  of  MCVs,  Bartels  and  Maddox  (1991) 
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concluded  that  “the  rapid  mesovortex  generation  observed 
can  be  explained  by  the  stretching  term  of  the  vorticity  equa¬ 
tion.”  This  conclusion  was  limited  to  those  MCVs  for  which 
a  visually  recognizable  signature  existed  in  satellite  imagery- 
preferentially  in  weak  shear.  Seven  years  later,  Weisman  and 
Davis  (1998)  argued  that  “for  systems  that  develop  in  an  en¬ 
vironment  with  weak-to-moderate  shear,  a  line-end  vortex 
pair  is  generated  primarily  via  the  tilting  of  horizontal  vor¬ 
ticity  generated  within  the  system-scale  cold  pool”  and  that 
“convergence  at  midlevels  enhances  Coriolis  rotation  over 
the  longer  term,  leading  to  the  preferred  development  of  a 
cyclonic  vortex”.  These  conclusions  certainly  apply  for  the 
storms  in  those  simulations.  However,  the  conclusions  may 
not  apply  to  all  storms  as  they  form  in  a  wide  variety  of 
background  conditions.  For  example,  storms  form  in  envi¬ 
ronments  with  vastly  different  thermodynamic  and  kinematic 
profiles,  and  they  also  evolve  in  response  to  atmospheric  mo¬ 
tions  over  a  broad  range  of  scales.  It  can  be  thought  that 
the  complex  array  of  possible  influencing  factors  in  storm 
development  would  broaden  the  spectrum  of  development 
characteristics  for  MCVs  beyond  the  traditional  “tilting  then 
stretching”  argument. 

This  wider  range  of  possibilities  applies  both  to  the  bud¬ 
gets  of  absolute  vertical  vorticity  and  to  the  respective  roles 
of  the  convective  and  stratiform  precipitation  regions  of 
MCSs.  Regarding  the  former,  for  example,  some  studies 
have  concluded  that  tilting  of  vorticity  is  the  primary  con¬ 
tributor  to  the  formation  of  MCVs  (Skamarock  et  ah,  1994; 
Weisman  and  Davis,  1998),  but  others  have  found  vortex  tube 
stretching  to  be  the  primary  vorticity  generation  mechanism 
at  mid-levels  (Brandes,  1990;  Bartels  and  Maddox,  1991; 
Johnson  and  Bartels,  1992;  Davis  and  Trier,  2002),  and  still 
others  have  suggested  that  both  mechanisms  may  be  impor¬ 
tant  at  different  stages  of  the  MCV  life  cycle  (Cram  et  ah, 
2002)  or  at  different  scales  of  atmospheric  motion  (Knievel 
and  Johnson,  2003).  Although  these  studies  may  appear  to 
be  in  conflict  with  each  other,  they  have  been  carried  out  in 
different  conditions  using  different  analysis  techniques  and 
different  model  domain  sizes  and  resolutions,  so  the  evolu¬ 
tion  need  not  be  the  same  in  all  cases.  Kirk  (2007)  made 
a  substantial  effort  to  resolve  the  apparent  confusion  among 
all  these  studies  by  using  a  phase  plot  method  to  track  the 
evolution  of  vertical  vorticity  budgets  in  different  MCVs.  To 
add  to  the  collection  of  cases  that  have  been  observed  and 
simulated,  and  with  the  intent  of  providing  some  clarifica¬ 
tion  of  the  MCV  formation  process  as  it  may  occur  in  moist 
baroclinic,  mid-latitude  environments,  we  perform  here  de¬ 
tailed  analyses  of  the  circulation  budget  of  two  simulated 
MCVs  over  their  lifetime,  from  genesis  to  maturity.  Addi¬ 
tionally,  our  idealized  simulations  will  start  with  a  basic  state 
in  thermal  wind  balance  and  include  the  effects  of  larger  scale 
baroclinic  environments,  something  which  was  not  included 
in  the  Davis  and  Weisman  (1994),  Skamarock  et  al.  (1994), 
Weisman  and  Davis  (1998),  or  Cram  et  al.  (2002)  simula¬ 
tions. 


As  a  part  of  this  goal,  we  will  provide  a  qualitative 
overview  of  the  PV  development  process  and  the  aggrega¬ 
tion  process  of  PV  elements  in  order  to  make  a  cursory  in¬ 
vestigation  into  the  role  of  PV  concentration  on  the  storm 
scale  (small,  intense,  convective  updrafts)  as  well  as  on  the 
MCV  scale.  Such  a  look  is  intended  to  qualitatively  illus¬ 
trate  the  processes  behind  the  more  quantitative  budgets  that 
appear  earlier  in  the  paper,  highlighting  the  role  of  vortic¬ 
ity  on  smaller  scales,  which  is  primarily  generated  by  tilting 
mechanisms,  on  the  development  of  vorticity  at  the  larger, 
MCV  scale  as  well  as  the  flux  convergence  of  planetary  vor¬ 
ticity  and  relative  vorticity  on  the  system-wide  scale.  Given 
the  linear  hodographs  employed  in  all  basic  states  used  to 
date  in  idealized  studies  of  MCV  dynamics,  one  would  ex¬ 
pect  tilting  to  generate  positive  and  negative  relative  vortic¬ 
ity  centers  in  nearly  equal  amounts  (e.g.,  Rotunno,  1981). 
When  tilting  occurs  in  an  environment  with  horizontal  vor¬ 
ticity,  due  to  a  convective  updraft,  centers  of  positive  and 
negative  vorticity  are  produced.  When  vertical  gradients  of 
diabatic  heating  act  upon  these  vorticity  centers,  they  pro¬ 
duce  a  concentration  of  PV  substance  (e.g.  Tory  et  al.,  2007) 
in  mid-levels  and  a  dilution  of  PV  substance  above,  resulting 
in  an  approximate  vertical  stacking  of  positive  and  negative 
PV  anomalies,  respectively.  The  evolution  of  these  positive 
and  negative  PV  centers  can  be  monitored  in  numerical  sim¬ 
ulations  to  discover  their  relative  contributions  to  the  MCV 
PV. 

The  magnitude  of  convective  versus  stratiform  processes 
contributing  to  the  PV  concentration  in  MCVs  is  less  well 
understood.  Many  studies  attribute  the  PV  concentration  to 
diabatic  heating  gradients  that  are  present  in  the  stratiform 
precipitation  region  of  an  MCS  (Hertenstein  and  Schubert, 
1991;  Johnson  and  Bartels,  1992;  Fritsch  et  al.,  1994;  Davis 
and  Trier,  2002).  This  traditional  view  is  supported  by  the 
fact  that  MCVs  are  consistently  observed  to  develop  in  the 
stratiform  portion  of  MCSs.  However,  the  diabatic  heating 
profiles  presented  by  Houze  (1997,  2004)  imply  that  the  PV 
concentration  may  come  from  either  convective  or  stratiform 
processes.  A  look  at  the  conceptual  model  of  the  MCS,  such 
as  presented  in  Fig.  14  of  Houze  (2004),  with  air  ascending 
in  convective  elements  at  the  leading  edge  of  the  MCS  then 
moving  rearward  to  the  stratiform  precipitation  region,  re¬ 
veals  that  the  PV  concentration  in  stratiform  and  convective 
areas  of  the  MCS  are  part  of  the  same  process  and  that  per¬ 
haps  the  PV  concentration  really  need  not  be  separated  into 
stratiform  and  convective  components.  Nevertheless,  MCSs 
do  have  varying  amounts  of  convective  and  stratiform  pre¬ 
cipitation  within  them,  and  the  purpose  of  our  analysis  is  to 
provide  evidence  that  convective  processes  may  be  a  more 
substantial  contributor  to  PV  concentration  than  has  been  tra¬ 
ditionally  thought. 

Deep  convection  often  occupies  a  relatively  small  area 
(by  percentage)  of  an  MCS,  but  the  diabatic  heating  rates 
and  vertical  mass  flux  per  unit  area  are  at  least  an  order  of 
magnitude  greater,  locally,  in  convective  than  in  stratiform 
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Table  1.  Selected  parameters  for  simulations. 


Control 

CAPE 

Basic  State 

Tcent  (K) 

293 

303 

Shearcent  (m  s“  1  km- 1 ) 

1.5 

1.5 

Width#  z  (degrees) 

60 

60 

CAPEcent  (Jkg-1) 

0 

2247 

CAPEniax  (Jkg-1) 

0 

4200 

Inner  domains 

Start  time  (hours) 

120 

120 

Move  frequency  (hours) 

D3 

- 

72 

D4 

24 

8 

Grid  dimensions 

D2 

94x76 

154x94 

D3 

229x172 

289x199 

D4 

271x211 

289x253 

Initial  vortex  radial  profile 

Guassian 

Isolated 

precipitation.  Additionally,  the  early  stages  of  the  life  cycle 
of  many  MCSs  tend  to  be  dominated  by  convection,  lead¬ 
ing  to  the  distinct  likelihood  that,  although  MCVs  are  rarely 
(if  ever)  observed  in  the  convective  region  of  an  MCS,  con¬ 
vection  can  be  the  primary  contributor  to  the  development  of 
PV  anomalies  that  are  eventually  revealed  in  the  stratiform 
precipitation  area.  Again,  since  most  of  the  mass  flux  is  pro¬ 
cessed  through  both  the  leading  convective  region  and  the 
stratiform  precipitation,  it  can  be  argued  that  both  may  be 
similar  contributors.  The  present  study  will  clarify  this  issue 
by  performing  a  detailed  analysis  of  the  diabatic  heating  pro¬ 
files.  These  diabatic  heating  profiles  should  reveal  the  quan¬ 
tity  of  mass  processed  through  the  convective  system  since 
diabatic  heating  rate  is  proportional  to  the  mass  flux. 

2  Simulations 

The  simulations  were  part  of  an  investigation  into  the  mech¬ 
anisms  for  the  development  of  MC  V s  in  weak  but  non-trivial 
baroclinic  environments  by  Conzemius  et  al.  (2007)  and  are 
further  described  in  that  study.  The  goal  of  the  original  study 
was  to  conduct  idealized  simulations  capable  of  coarsely  re¬ 
solving  deep  convective  processes  in  an  environment  with 
relatively  weak  shear,  whereby  the  dynamics  of  MCV  main¬ 
tenance  and  growth  could  be  examined,  but  still  encompass¬ 
ing  an  area  large  enough  that  permits  larger  scale  baroclinic 
development  as  part  of  the  MCV  lifecycle. 

For  this  study,  we  conducted  two  simulations,  the  pri¬ 
mary  difference  between  which  was  the  value  of  ambient 
CAPE.  The  first  simulation  was  designated  the  control  sim¬ 
ulation  and  started  with  a  basic  state  thermodynamic  profile 
that  was  neutral  to  the  moist  ascent  of  boundary  layer  air 
parcels.  By  definition,  this  neutral  atmosphere  had  no  CAPE, 
but  CAPE  developed  during  the  simulation,  reaching  maxi¬ 
mum  values  of  around  700  J  kg-1.  The  second  simulation 
had  2000  J  kg-1  CAPE  in  the  center  of  the  baroclinic  zone 
and  was  simply  designated  as  the  CAPE  simulation. 


We  performed  both  simulations  using  the  MM5  model 
(Grell  et  al.,  1994),  which  is  a  limited  area,  non-hydrostatic, 
sigma  coordinate  model.  We  performed  the  simulations  on 
an  1 1  000  km  by  5000  km  domain  using  an  outer,  coarse- 
resolution  grid  with  a  90  km  grid  interval  and  three  nested 
inner  grids  whose  intervals  were  30,  10,  and  3.3  km.  The  ini¬ 
tial  state  consisted  of  an  Eady-type  baroclinic  background, 
in  which  all  variables  are  a  function  of  latitude  and  pressure 
only.  The  initial  state  differs  from  the  Eady  basic  state,  how¬ 
ever,  in  that  a  stratosphere  was  included,  and  the  shear  was 
half  the  value  of  the  Eady  basic  state.  The  lateral  boundary 
conditions  were  prescribed  to  be  the  initial,  basic  state  con¬ 
ditions,  the  upper  boundary  condition  was  a  gravity  wave  ra¬ 
diation  condition,  and  the  bottom  boundary  conditions  were 
no  slip  but  free  of  sensible  or  latent  heat  fluxes.  Addition¬ 
ally,  there  was  no  infrared  or  shortwave  radiation.  Without 
heat  fluxes  or  radiation,  there  was  no  diurnal  cycle  in  the 
experiments.  Further  details  and  reasoning  behind  these  ex¬ 
periments  are  provided  in  Conzemius  et  al.  (2007).  The  basic 
parameters  of  the  simulations  are  listed  in  Table  1 . 

In  order  to  initiate  an  MCS,  we  inserted  a  warm  core, 
low-level  vortex,  with  zero  PV  perturbation  in  its  interior  (to 
avoid  preconditioning  the  simulation  with  interior  PV)  into 
the  background  state.  The  vortex  had  a  radial  profile  of  po¬ 
tential  temperature  perturbation  with  maximum  magnitude 
of  3.0 K  at  the  surface  and  decayed  with  height.  The  radial 
profile  was  in  the  shape  of  a  Gaussian  function  for  the  con¬ 
trol  simulation  and  in  the  shape  of  the  second  derivative  of 
the  Gaussian  profile  (i.e.  “isolated”)  for  the  CAPE  simula¬ 
tion.  As  it  turns  out,  the  primary  effect  of  the  initial  vortex 
in  the  simulations,  as  discussed  in  Conzemius  et  al.  (2007), 
is  to  excite  a  spectrum  of  moist  baroclinic  waves.  The  MCS 
develops  just  to  the  east  of  the  surface  low  that  is  associ¬ 
ated  with  the  most  rapidly  amplifying  baroclinic  wave.  The 
isolated  profile  was  chosen  for  the  CAPE  simulation  in  an 
attempt  to  minimize  the  impact  of  the  moist  baroclinic  in¬ 
stability.  However,  since  the  MCS  development  inevitably 
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was  tied  to  the  baroclinic  processes,  this  change  served  only 
to  slow  the  large  scale  moist  baroclinic  growth  and  thereby 
delay  the  MCS  development. 

Conzemius  et  al.  (2007)  describe  the  evolution  of  these 
simulations  in  greater  detail,  so  we  will  only  summarize  here. 
In  the  control  simulation,  an  MCS  developed  approximately 
126  h  into  the  simulation,  within  an  intensifying,  moist  baro¬ 
clinic  cyclone.  An  MCV  developed  within  the  MCS,  and 
the  MCS  precipitation  gradually  became  more  stratiform  in 
nature  as  the  simulation  proceeded.  In  the  CAPE  simula¬ 
tion,  the  MCS  developed  approximately  150  h  into  the  simu¬ 
lation  (this  was  delayed  due  to  the  choice  of  an  isolated  initial 
vortex  temperature  radial  profile),  but  when  the  growth  oc¬ 
curred,  it  was  much  more  rapid  and  extensive.  The  MCV  that 
developed  was  much  stronger  as  well,  with  the  precipitation 
shield  remaining  convective  much  longer  into  the  simulation. 


3  Diagnostic  calculations 
3.1  Vorticity  budget  calculations 


Similar  to  the  procedure  expounded  in  Haynes  and  McIntyre 
(1987)  and  the  methodology  outlined  in  Weisman  and  Davis 
(1998),  Cram  et  al.  (2002)  and  Davis  and  Trier  (2002),  we 
use  the  flux  form  of  the  vorticity  equation  and  calculate  area- 
and  temporal-averages  of  the  circulation  tendency: 
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where  / a  is  the  absolute  vorticity  defined  as 
C,  =  /  +  dv/dx  —  du/dy,  and  F  and  G  represent  the 
horizontal  components  of  the  effects  of  additional  forces, 
such  as  friction,  which  are  neglected  in  this  analysis  for 
simplicity.  The  neglect  of  frictional  forces  can  be  justified  by 
the  fact  that  the  initial  near-surface  winds  are  very  weak  and 
is  further  supported  by  the  comparisons  between  vorticity 
budget  profiles  and  actual  vorticity  change  presented  in 
Sect.  4.1. 

In  addition  to  avoiding  the  inherent  problem  of  measuring 
the  small  difference  between  large  terms  that  can  occur  when 
using  the  material  form  of  the  vorticity  equation  (Haynes  and 
McIntyre  1987),  the  strength  of  this  method  lies  in  the  fact 
that,  through  Gauss’s  theorem,  the  tendency  of  the  circula¬ 
tion  for  any  enclosed  area  fixed  in  space  can  be  written  in 
terms  of  the  line  integral  of  the  flux  component  normal  to  the 
boundary  of  the  area  (Davis  and  Trier,  2002): 


where  C  is  the  absolute  circulation,  A  represents  the  area 
encompassed  by  the  closed  loop,  and  u  the  horizontal  ve¬ 
locity.  We  will  refer  to  the  first  term  under  the  integral  as 
the  horizontal  convergence  of  absolute  vorticity  (horizontal 


flux  term).  A  comparison  between  the  flux  and  the  mate¬ 
rial  forms  of  the  vorticity  equation  reveals  that  the  horizontal 
flux  term  in  Eq.  (2)  is  equivalent  to  the  horizontal  advection 
term  plus  the  stretching  term  in  the  material  form  of  the  vor¬ 
ticity  equation.  The  second  term  is  an  amalgamation  of  the 
vertical  advection  of  vertical  vorticity  and  the  tilting  of  hor¬ 
izontal  vorticity  into  vertical  vorticity  and  will  be  hereafter 
denoted  simply  as  the  tilting-like  term  (see  Tory  et  al.  (2007) 
for  clarification  of  this  terminology).  For  comparison,  Davis 
and  Trier  (2002)  refer  to  these  terms  as  stretching  and  tilting, 
respectively. 

The  corresponding  form  of  the  PV  equation  is  most  conve¬ 
niently  written  in  theta  (potential  temperature)  coordinates: 

A 

where  Q  is  PV  (m2  s-1  Kkg-1),  cr  the  isentropic  density 
(m-2  K_  1  kg),  and  6  the  material  rate  of  change  of  poten¬ 
tial  temperature.  In  theta  coordinates,  the  PV  substance  crQ 
is  simply  dv/dx  —  du/dy  +  /,  so  one  can  see  a  striking  sim¬ 
ilarity  between  Eqs.  (2)  and  (3)  (the  zonal  and  meridional 
partial  derivatives  of  course  being  taken  on  surfaces  of  con¬ 
stant  pressure  and  potential  temperature,  respectively).  Tory 
et  al.  (2007)  and  Haynes  and  McIntyre  (1987)  discuss  at 
length  the  corresponding  PV  equation  in  differential  form. 
The  divergence  of  the  second  (non-advective)  flux  term  can 
be  split  into  two  terms,  one  of  which  describes  tilting-like 
effects  and  the  other  describing  vertical  advection-like  ef¬ 
fects.  The  tilting-like  and  vertical  advection-like  terms  tend 
to  largely  balance  each  other,  typically  leaving  the  first  term 
in  the  line  integral  as  the  dominant  term  on  the  right  hand 
side  of  Eq.  (3).  Thus,  we  would  expect  an  increase  in  PV 
substance  over  any  area  to  be  mostly  due  to  flux  convergence 
(i.e.  a  concentration)  of  PV.  This  concentration  of  PV  sub¬ 
stance  occurs  mostly  at  low  levels,  driven  by  diabatic  heating 
within  the  MCS.  At  upper  levels,  we  would  expect  a  nega¬ 
tive  PV  anomaly  to  form  as  flux  divergence  dilutes  the  PV 
substance. 

3.2  Areas  and  times  of  analysis  on  innermost  grid 

The  rectangle  forming  the  contour  of  integration  was  chosen 
to  be  as  close  to  the  vortex  center  while  encompassing  all  sig¬ 
nificant  cyclonic  relative  vorticity.  In  practice,  the  southern 
line  segment  was  placed  along  the  maximum  850-700  hPa 
rear  inflow  into  the  convective  system,  while  the  western  and 
northern  line  segments  were  placed  as  close  as  possible  to 
the  cyclonic  vorticity  maximum  but  in  regions  encountering 
little  or  no  convection.  To  allow  for  the  movement  of  the  con¬ 
vective  system  within  the  averaging  time  interval,  we  placed 
the  eastern  line  segment  far  enough  ahead  of  the  convective 
line  that  the  precipitation  did  not  advance  beyond  the  eastern 
line  until  after  the  end  of  the  averaging  time  interval. 

Four  time  periods,  each  of  288  min  duration,  were  chosen 
for  the  vorticity  budget  analysis  of  the  control  simulation. 
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Table  2.  Analysis  box  characteristics. 


Simulation 

Time  interval  (h) 

Dimensions  (grid  cells) 

Area  (km2) 

Control 

A 

124.8-129.6 

65x45 

32500 

B 

129.6-134.4 

110x60 

73  300 

C 

134.4-139.2 

140x70 

109  000 

D 

170.4-175.2 

175x110 

214000 

CAPE 

A 

152-156 

100x56 

62200 

B 

156-160 

140x110 

171000 

C 

160-164 

145x130 

209  000 

D 

168-172 

230x160 

409  000 

«f  Bl  71  BI  ST 

4  « 

3  V  8  W  IV 

71  av  SI  4  V 

31  8  V  I  W  0  IK 

Fig.  1.  Model-derived  reflectivity  (dBz;  shading),  horizontal  flow  (barbs;  knots)  at  700  hPa,  and  analysis  boxes  at  the  following  times  in  the 
control  simulation:  (a)  129.6  h;  (b)  132.0  h;  (c)  136.8  h;  and  (d)  172.8  h.  The  equivalent  radar  reflectivity  factor  for  rain,  Ze  (mm6  m-3),  is 
computed  following  Fovell  and  Ogura  (1988)  as  Ze  =  120akNok~^  ( P/Pw )2,  where  a  is  the  nondimensional  ratio  of  backscattering  for  the 
reflecting  substance  (1  for  water,  0.213  for  ice),  k=  1018  is  the  conversion  factor  from  m6  to  the  units  of  Ze,  Na  is  the  intercept  parameter 
of  the  raindrop  size  distribution,  p  the  density  of  the  substance  of  interest  (water  or  ice),  and  pw  the  density  of  water,  k  is  the  slope  of  the 
raindrop  size  distribution:  k  =  (npwN0l  paqra )3'  ,  where  pa  is  the  density  of  dry  air,  and  qra  the  rainwater  mixing  ratio.  The  reflectivity  is 
expressed  in  decibels,  or  dBZ,  where  dBZ  =  lOlogjo(Ze). 
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Fig.  2.  Model-derived  reflectivity  (dBz;  shading),  horizontal  flow  (barbs;  knots)  at  700  hPa,  and  analysis  boxes  at  the  following  times  in  the 
CAPE  simulation:  (a)  156  h;  (b)  158  h;  (c)  162  h;  and  (d)  170  h. 


The  chosen  time  periods  encompass  1 .  initiation  and  organi¬ 
zation  of  deep  convection  accompanied  by  intensification  of 
the  baroclinic  cyclone,  2.  expansion  of  the  stratiform  precip¬ 
itation  region  and  slight  weakening  of  convection  (very  little 
deepening  occurs  during  this  phase),  3.  a  period  of  essen¬ 
tially  no  deepening  of  the  baroclinic  cyclone,  and  4.  a  later 
period,  characterized  by  a  resumption  of  deepening.  The  line 
integral  was  calculated  at  each  model  level  on  the  innermost 
3.3-km  domain  at  8  min  intervals,  and  time  averaging  was 
performed  over  the  288-min  analysis  period. 

For  the  CAPE  simulation,  the  time  averaging  periods  were 
240  min  each.  The  first  period  corresponds  to  the  develop¬ 
ment  and  organization  of  the  MCS  and  intensification  of  the 
baroclinic  cyclone.  The  next  two  periods  correspond  to  the 
continued  expansion  of  the  stratiform  precipitation  area,  dur¬ 
ing  which  time  the  baroclinic  cyclone  did  not  deepen  (at  least 
in  terms  of  surface  pressure)  significantly.  The  final  period 


saw  a  resumption  of  rapid  deepening  of  the  baroclinic  cy¬ 
clone. 

The  placement  of  the  analysis  boxes  is  shown  for  the  con¬ 
trol  simulation  in  Fig.  1  and  for  the  CAPE  simulation  in 
Fig.  2.  In  Table  2,  we  provide  information  about  the  char¬ 
acteristics  of  the  analysis  boxes.  Because  of  the  increasing 
size  of  the  circulating  flow  in  the  simulations,  we  increased 
the  size  of  the  analysis  box  between  successive  analysis  in¬ 
tervals  in  order  to  encompass  the  cyclonic  relative  vorticity. 

3.3  Convective  versus  stratiform  partitioning  method 

By  partitioning  the  updraft  mass  flux  into  a  convective  and 
a  stratiform  portion,  one  can  assess  the  contributions  of  con¬ 
vective  and  stratiform  areas  to  the  concentration  of  absolute 
vorticity  or  PV  within  the  system.  The  convective-stratiform 
partitioning  method  we  used  is  similar  to  that  used  in  Tao 
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et  al.  (1993)  and  Braun  (2009)  except  that  the  thresholds  for 
defining  convective  areas  are  somewhat  higher.  The  parti¬ 
tioning  was  done  on  a  two-dimensional,  grid  column  basis. 
Grid  columns  were  designated  convective  if  the  surface  pre¬ 
cipitation  rate  was  twice  the  average  precipitation  rate  of  the 
immediate,  surrounding  four  cells.  In  such  cases,  the  sur¬ 
rounding  four  cells  were  also  designated  as  convective.  Sec¬ 
ondly,  all  grid  cells  with  a  rainfall  rate  exceeding  25  mmh-1 
were  assigned  to  the  convective  bin.  In  some  situations, 
strong  convective  updrafts  may  not  have  rain  reaching  the 
surface,  so  if  vertical  velocity  at  any  grid  point  in  the  col¬ 
umn  exceeded  5  ms-1,  or  if  the  cloud  liquid  water  content 
exceeded  0.5  g  kg-1  below  the  freezing  level,  then  the  col¬ 
umn  was  designated  as  convective.  For  the  remaining  (non- 
convective)  grid  columns,  stratiform  areas  were  defined  as 
those  with  a  precipitation  rate  of  at  least  0.1  mmh-1.  Grid 
cells  within  the  analysis  box  meeting  none  of  these  criteria 
were  assigned  to  the  non-precipitating  category.  The  areas 
and  time  periods  of  analysis  were  the  same  as  for  the  vortic¬ 
ity  budget  analysis. 

Since  the  tendency  of  area  average  potential  vorticity  sub¬ 
stance  is  driven  by  the  horizontal  flux  convergence  of  PV 
substance  (Eq.  3),  and  because  the  absolute  vorticity  ten¬ 
dency  is  also  affected  by  horizontal  flux  convergence  (Eq.  2), 
one  can  invoke  mass  conservation  to  relate  the  flux  conver¬ 
gence  of  either  quantity  to  the  updraft  mass  flux  within  the 
MCS.  Since  the  vertical  mass  flux  is  directly  related  to  the 
diabatic  heating  rate  in  the  system  (we  assume  that  the  dia- 
batic  heating  rate  is  essentially  balanced  by  the  upward  trans¬ 
port  of  smaller  potential  temperature),  it  is  equally  useful  to 
use  diabatic  heating  rate  as  a  proxy  for  mass  flux.  Our  ini¬ 
tial  calculations  demonstrated  that  the  mass  flux  profiles  (not 
shown)  were  qualitatively  very  similar  to  the  profiles  of  dia¬ 
batic  heating  rate.  Similar  to  Kirk  (2007),  the  diabatic  heat¬ 
ing  rate  was  calculated  by  taking  the  residual  of  the  thermo¬ 
dynamic  equation  (time  tendency  minus  advection)  and  av¬ 
eraging  over  horizontal  planes  within  the  analysis  area.  We 
assumed  diffusion  effects  on  the  thermodynamic  equation  to 
be  negligible. 

4  Sources  of  absolute  vorticity  and  PV 
4.1  Vorticity  budget  profiles 

During  most  of  the  time  periods  of  our  analysis,  the  tilting- 
like  term  is  a  contributor  to  the  creation  of  positive  vertical 
vorticity  in  the  lower  troposphere  up  to  about  700  hPa,  but 
its  calculated  contribution  is  sensitive  to  the  placement  of  the 
rectangle  used  to  calculate  the  loop  integrals.  The  tilting-like 
term  generally  becomes  large  during  brief  intervals  when  the 
southern  line  segment  intersects  individual  convective  cells, 
which  make  strong  contributions  to  the  vorticity  budget  on 
scales  smaller  than  the  size  of  the  MCS  itself  (Knievel  and 
Johnson  2003).  The  four  hour  integration  time  smoothes 


out  these  contributions,  however.  Analyses  on  the  coarser 
grids  (not  shown),  which  do  not  resolve  individual  convec¬ 
tive  cells,  reveal  similar-looking  profiles  of  vorticity  budget 
terms,  suggesting  some  MCS  scale  tilting  as  well.  Never¬ 
theless,  the  circulation  budget  during  the  earliest  analysis 
period  (Fig.  3a)  is  dominated  by  the  horizontal  flux  conver¬ 
gence  of  absolute  vorticity.  This  finding  refutes  the  “tilting 
then  stretching”  argument  proposed  by  Weisman  and  Davis 
(1998)  for  this  particular  case.  As  discussed  in  the  intro¬ 
duction,  the  initial  state  differs  from  the  Weisman  and  Davis 
(1998)  simulations  due  to  its  zero  initial  CAPE  and  its  inclu¬ 
sion  of  baroclinic  processes. 

This  finding  is  also  somewhat  different  from  that  of  Kirk 
(2007),  where  the  tilting-like  term  (vertical  momentum  flux 
curl)  contributed  most  strongly  in  the  mid-troposphere.  The 
analysis  technique  we  employed  here,  however,  is  differ¬ 
ent  from  that  used  in  Kirk  (2007),  and  consequently,  it  is 
rather  difficult  to  attribute  the  same  meaning  to  the  results. 
In  particular,  our  analysis  box  was  fixed  in  time,  whereas  in 
Kirk  (2007),  the  box  moves  with  the  MCV.  Additionally,  the 
Kirk  analysis  box  is  72kmx72km  (an  area  of  5184km2), 
whereas  our  analysis  boxes  ranged  in  size  from  32  500  km2 
to  409  000  km2. 

On  the  scale  of  the  MCS,  the  horizontal  flux  convergence 
is,  by  far,  the  dominant  mechanism  responsible  for  the  spin 
up  of  mid-level  cyclonic  absolute  vorticity.  The  horizontal 
flux  convergence  in  Fig.  3a  and  b  is  strongest  at  low  lev¬ 
els,  consistent  with  a  divergence  profile  associated  with  deep 
convection  (Houze  1997)  and  also  consistent  with  the  find¬ 
ings  of  Kirk  (2007). 

By  t  =  129.6  h  (Fig.  3b),  a  closed  circulation  is  evident  at 
700  hPa  (not  shown).  Between  t  =  124.8  and  t  =  134.4h  into 
the  simulation,  the  mid-level  circulation  strengthens  substan¬ 
tially.  This  pattern  is  consistent  with  the  re-intensification 
of  the  27-28  May  1998  MCV  investigated  by  Davis  and 
Trier  (2002).  The  stratiform  precipitation  is  also  expanding 
rapidly,  and  during  this  time,  the  horizontal  flux  convergence 
of  absolute  vertical  vorticity  is  reflective  of  the  stratiform  di¬ 
vergence  profile  suggested  by  Houze  (1997),  with  its  flux 
convergence  maximum  at  higher  altitude  than  in  Fig.  3a.  In 
the  ensuing  time  period  (Fig.  3c),  the  circulation  tendency  is 
strongly  positive  in  the  upper  troposphere,  due  to  the  hori¬ 
zontal  flux  divergence  of  negative  relative  vorticity. 

At  a  much  later  time  (Fig.  3d),  when  the  surface  low  un¬ 
dergoes  a  steadier  strengthening,  the  circulation  budget  is  al¬ 
most  entirely  dominated  by  the  horizontal  flux  term.  The 
dominance  of  this  term  shows  that  the  mean  secondary  circu¬ 
lation  associated  with  the  vortex  (convergence  at  low  levels, 
rising  in  the  vicinity  of  the  vortex,  and  divergence  at  upper 
levels)  drives  the  strengthening  of  the  circulation,  but  this  re¬ 
sult  might  also  be  an  artifact  of  the  analysis  area.  Due  to  the 
size  and  orientation  of  the  MCS,  it  was  difficult  to  place  the 
southern  portion  of  the  line  integral  along  the  maximum  rear 
inflow  into  the  system.  The  system  is  also  somewhat  less 
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(a)  /=  124.8  h  to/=  129.6  h  (b)/=  129.6  h  to/=  134.4  h 


(c)  /  =  134.4  h  to  t  =  1 39.2  h  (d)  /  =  1 711.4  >■  to  /  =  1 75.2  h 


Fig.  3.  Circulation  budget  analysis  expressed  as  area  average  vorticity  change  for  a  rectangle  surrounding  the  vortex  in  the  control  simulation 
on  the  3. 3 -km  grid  for  the  time  periods:  (a)  124.8h  to  129.6 h;  (b)  129.6h  to  134.4h;  (c)  134.4h  to  139.2 h;  and  (d)  170.4h  to  175.2 h, 
respectively.  The  dotted  line  indicates  the  line  integral  of  horizontal  flux  convergence  of  absolute  vorticity;  the  dashed  line  is  the  correspond¬ 
ing,  tilting-like  term;  the  solid  black  line  is  the  total  of  the  horizontal  and  tilting-like  terms;  and  the  solid  gray  line  is  the  actual  change  over 
the  period  of  integration. 


convective  at  this  time.  Section  4.2  will  explore  the  convec¬ 
tive  and  stratiform  partitioning  in  greater  detail. 

The  profiles  of  the  vorticity  budget  terms  at  early  times 
(Fig.  3a,  b)  resemble  those  in  Davis  and  Galarneau  (2009), 
particularly  with  regard  to  the  peak  of  flux  convergence  in 
the  middle  troposphere  and  the  dominance  of  flux  conver¬ 
gence  in  general  (N.b.  our  flux  convergence  term  is  the  sum 
of  the  eddy  flux  and  stretching  terms  defined  in  Davis  and 
Galarneau).  In  both  analyses,  the  budget  terms  were  calcu¬ 
lated  using  rectangles  of  similar  size,  the  MCVs  occurred 
within  similar  values  of  deep  layer  shear,  and  baroclinic  ef¬ 
fects  were  included  in  the  simulations.  One  particular  differ¬ 
ence  in  the  vorticity  budget  terms,  however,  lies  in  the  tilting- 
like  term  in  the  lower  troposphere,  which  is  negative  in  the 
Davis  and  Galarneau  analyses  but  positive  in  ours.  This  dif¬ 
ference  may  be  due  to  the  eastern  and  northern  portions  of 
the  box  being  located  in  precipitation  areas  in  the  Davis  and 


Galarneau  cases.  Further  analysis  is  required  to  uncover  the 
root  of  these  differences,  but  the  overall  result  is  the  same. 
That  is,  tilting  is  secondary  to  stretching  in  the  production  of 
absolute  vorticity. 

In  the  CAPE  simulation  (Fig.  4),  the  evolution  of  the  vor¬ 
ticity  budget  is  somewhat  different.  In  the  initial  time  period 
(Fig.  4a),  the  horizontal  flux  convergence  term  is  responsi¬ 
ble  for  all  the  spin-up  of  absolute  vorticity  below  the  500  hPa 
level,  except  for  the  lowest  100  hPa  (where  cold  pool  effects 
occur  -  see  below).  Once  again,  the  overall  profile  and  domi¬ 
nance  of  the  flux  convergence  term  is  similar  to  the  results  of 
Davis  and  Galarneau  (2009)  during  this  initial  analysis  pe¬ 
riod.  Thereafter  (Fig.  4b  and  c),  the  flux  convergence  term 
actually  becomes  a  negative  contributor  to  the  absolute  vor¬ 
ticity  below  about  600  hPa,  and  the  spin  up  of  absolute  vortic¬ 
ity  is  provided  solely  by  the  tilting-like  term.  It  appears  that 
the  cold  pool  generation  is  influencing  the  flux  convergence 
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(a) /=  152h  lo/=  156h  (b) /  =  156h  to/=  I60h 


(c)  t  =  1 60h  lo  /  =  I  (>4h 


(d)  I  -  I68h  to  /=  1 721i 


Fig.  4.  Circulation  budget  analysis  expressed  as  area  average  vorticity  change  for  a  rectangle  surrounding  the  vortex  on  the  3.3-km  grid  in 
the  CAPE  simulation.  The  time  periods  are  as  follows:  (a)  152h  to  156 h;  (b)  156h  to  160 h;  (c)  160h  to  164 h;  and  (d)  168h  to  172h, 
respectively.  For  notation,  see  Fig.  3. 


during  this  time  as  the  plots  of  lOOOhPa  temperature  show  a 
significant  cold  anomaly  (as  much  as  10°C)  spreading  radi¬ 
ally  outward.  In  mid  levels,  however,  the  flux  convergence  is 
a  contributor  to  the  intensification  of  cyclonic  absolute  vor¬ 
ticity.  At  the  very  end  of  the  time  period  of  analysis  (Fig.  4d), 
the  horizontal  flux  convergence  again  dominates  the  abso¬ 
lute  vorticity  budget  (again,  similar  to  Davis  and  Galarneau, 
2009).  By  this  time,  however,  the  area  of  integration  is  so 
large  that  any  mesoscale  tilting  contributes  less  to  the  area- 
integrated  vorticity. 

4.2  Stratiform  and  convective  contributions  to  vortex 
intensification 

A  time  series  of  the  convective  and  stratiform  areas,  within 
the  analysis  box,  during  the  analysis  intervals  in  the  control 
simulation  (not  shown)  reveals  that  the  convection,  as  defined 
in  Sect.  3.3,  never  occupies  more  than  ten  percent  of  the  area 
of  analysis.  Convection  initially  comprises  100  percent  of 


the  total  precipitation  area  (defined  here  as  the  sum  of  con¬ 
vective  and  stratiform  areas),  but  as  the  stratiform  precipita¬ 
tion  area  expands,  the  convective  area  decreases  to  as  little  as 
7  percent  of  the  precipitation  area  during  the  second  analysis 
period.  For  the  remaining  analysis  periods,  the  convection 
recovers  to  between  20  and  30  percent  of  the  total  precipita¬ 
tion  area.  By  the  end  of  the  final  analysis  period  (nearly  two 
days  after  the  beginning  of  the  first),  the  total  precipitation 
area  has  reached  60  000  km2. 

In  the  CAPE  simulation,  the  precipitation  area  is  consider¬ 
ably  larger,  reaching  a  maximum  of  177  000  km2.  The  con¬ 
vection  starts  out  occupying  89  percent  of  the  total  precipita¬ 
tion  area,  but  by  r  =  155.5  h,  it  drops  to  less  than  20  percent 
and  then  remains  between  10  and  20  percent  for  the  remain¬ 
ing  analysis  times.  The  convective  area  expands  throughout 
the  simulation  and  reaches  over  20000  km2  by  about  19  h 
after  the  start  of  the  first  analysis  period.  The  stratiform 
area  undergoes  a  similar  expansion  throughout  the  simula¬ 
tion.  In  the  CAPE  simulation,  the  convection  provides  a 
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(a)  t=  124.8  -  129.6  h  (d)/=  152  -  156  li 


-8  -4  0  4  8  -80  -40  0  40  80  120 


Diabatic  Heating  (K  day-1)  Diabalic  Heating  (K.  day') 


(b) /=  154.4  -  139.2  h  (e)  /=  160  -  164  h 


-8  -4  0  4  8  -40  -20  0  20  40 

Diabatic  Heating  (K  day')  Diabatic  Heating  (K.  day"1) 


(C)  /  =  170.4-  175.2  h  (f)  1=108- 172  h 


-8  -4  0  4  8  -20  -10  0  10  20  30 

Diabatic  Heating  (K.  day')  Diabatic  Heating  (K.  day') 


Fig.  5.  Area-averaged  diabatic  heating,  showing  contributions  from  convective  and  stratiform  precipitation  regions  in  the  control  simulation, 
for  a  box  surrounding  the  vortex  on  the  3.3-km  grid,  at  simulation  time:  (a)  124.8  h  to  129.6  h;  (b)  134.4  h  to  139.2  h;  and  (c)  170.  h  to 
175.2 h,  and  for  the  CAPE  simulation  at  (d)  152h  to  156 h;  (e)  160  h  to  164h;  and  (f)  168h  to  172 h.  The  area  average  heating  rate  pertains 
to  the  analysis  box,  which  is  centered  on  the  area  of  the  circulation  and  has  its  southern  edge  roughly  along  the  line  of  maximum  rear  inflow 
into  the  system. 
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much  greater  contribution  to  the  MCV  diabatic  heating  pro¬ 
file. 

Animations  of  relative  vorticity  (not  shown)  suggest  the 
dominance  of  tilting  at  smaller  scales  (e.g.  Knievel  and  John¬ 
son  2003),  but  on  an  MCV  scale,  the  horizontal  flux  con¬ 
vergence  is,  by  far,  the  dominant  contributor  to  the  vorticity 
change. 

The  evolution  of  the  diabatic  heating  profiles  in  the  con¬ 
trol  simulation  is  generally  reflective  of  the  evolution  of  the 
MCS  (see  Fig.  5).  In  the  initial  time  period  (Fig.  5a),  the  peak 
heating  is  nearly  equally  shared  by  convective  and  stratiform 
precipitation  regions.  Since  the  cooling  in  the  stratiform 
area  is  confined  to  the  lowest  150  hPa,  the  heating  gradient 
is  stronger  in  the  stratiform  region,  implying  the  stratiform 
precipitation  is  playing  the  larger  role  in  PV  concentration  in 
the  lower  troposphere.  By  r  =  1 34.4  h  (Fig.  5b),  the  stratiform 
precipitation  region  (the  non-convective  area  experiencing  a 
rainfall  rate  greater  than  0.1  mm  h-1)  has  expanded  to  oc¬ 
cupy  80  percent  of  the  precipitation  area.  Again,  the  diabatic 
cooling  in  the  lower  troposphere  provides  for  a  deeper  layer 
of  stratiform  PV  concentration.  At  this  point  in  the  simu¬ 
lation,  the  distinction  between  the  stratiform  and  convective 
profiles  is  consistent  with  the  characteristic  profiles  shown  in 
Fig.  3  of  Houze  (1997). 

During  the  final  time  period  (Fig.  5c),  as  the  system 
reaches  a  regime  of  steady  intensification,  the  diabatic  heat¬ 
ing  increases,  and  most  of  this  diabatic  heating  (and  there¬ 
fore  the  PV  concentration)  is  accomplished  by  the  stratiform 
precipitation  areas  within  the  system.  By  this  point,  the  di¬ 
abatic  heating  probably  reflects  the  moist  baroclinic  intensi¬ 
fication  at  least  as  much  as  it  is  a  sign  of  MCV  dynamics. 
Indeed,  Conzemius  et  al.  (2007)  show  that,  during  this  stage 
of  intensification,  the  conversion  of  mean  state  available  po¬ 
tential  energy  (APE)  to  eddy  APE  increases  in  tandem  with 
the  diabatic  production  of  eddy  APE.  Nevertheless,  the  MCS 
dynamics  are  a  contributor  to  this  process  by  providing  di¬ 
abatic  support  to  the  baroclinic  cyclone  intensification,  and 
balanced  MCV  dynamics  play  a  part  in  the  maintenance  of 
convection  (Conzemius  et  al.,  2007;  Trier  and  Davis,  2002). 
It  is  noteworthy  that,  during  this  final  period,  the  vertical  gra¬ 
dients  in  diabatic  heating  are  largest  in  the  lower  troposphere 
(much  like  the  convective  profile  is,  see  Houze,  2004),  yet  the 
precipitation  is  mostly  stratiform.  This  may  be  symptomatic 
of  the  fact  that  the  system  has  moved  north  into  areas  where 
the  mean  state  has  a  shallower  troposphere.  Additionally,  the 
system  has  grown  upscale  and  therefore  no  longer  resembles 
an  MCS. 

In  the  CAPE  simulation  (Fig.  5d  through  f),  the  diabatic 
heating  profile  is  markedly  different,  and  convective  heating 
plays  the  leading  role  in  PV  concentration.  At  first  (Fig.  5d), 
convection  contributes  essentially  all  of  the  diabatic  heat¬ 
ing,  and  the  diabatic  heating  gradient  is  strongest  between 
900  hPa  and  600  hPa.  The  stratiform  contribution  to  PV 
concentration  is  essentially  limited  to  a  weak  layer  between 
700  hPa  and  450  hPa.  As  the  stratiform  area  rapidly  expands 


(Fig.  5e),  its  contribution  to  the  total  becomes  more  signif¬ 
icant,  yet  with  a  diabatic  heating  difference  of  25  K  day-1 
over  the  700  hPa  to  500  hPa  layer,  it  still  produces  less  PV 
than  the  convective  region  (nearly  40  K  day- 1  difference  in 
heating  rate  between  600  hPa  and  the  surface).  It  is  thus  ev¬ 
ident  that  the  deep  convection  accomplishes  the  bulk  of  the 
flux  convergence  of  PV  substance  (see  Eq.  3)  until  the  very 
last  period  of  analysis  (Fig.  5f),  when  the  stratiform  diabatic 
heating  profile  becomes  more  consistent  with  the  idealized 
profile  shown  in  Houze  (1997)  and  is  more  favorable  for  con¬ 
tributing  to  PV  in  the  mid-troposphere.  In  this  final  analysis 
period,  the  convective  and  stratiform  heating  gradients  are 
nearly  equal. 

Overall,  convection  is  the  dominant  PV  concentration 
mechanism  in  the  CAPE  simulation.  But  clearly,  one  3.3  km 
simulation  cannot  prove  that  the  PV  concentration  occurs 
mostly  in  convective  updrafts  within  an  MCS.  In  fact,  the 

3.3  km  grid  interval  is  relatively  coarse  for  resolving  pro¬ 
cesses  essential  to  individual  convective  updrafts  (e.g.  Bryan 
et  al.,  2003).  However,  we  contend  that  the  simulation  per¬ 
forms  reasonably  well  at  resolving  processes  on  the  scale  of 
the  MCS,  and  the  results  are  believed  usefully  suggestive  of 
the  fact  that  much  of  the  PV  concentration  is  accomplished 
by  deep  convection.  While  the  deep  convective  dominance 
of  PV  concentration  is  very  evident  in  the  present  analysis, 
the  verification  of  these  results  also  requires  a  comprehensive 
sampling  of  the  observed  lifecycle  of  MCVs.  The  Bow  Echo 
and  MCV  Experiment  (BAMEX;  Davis  et  al.,  2004)  has  pro¬ 
vided  a  wealth  of  observational  datasets  that  can  be  used  to 
explore  this  problem  further. 

4.3  Trajectory  analysis  of  the  PV  evolution 

In  order  to  further  illustrate  the  evolution  of  PV  in  the  MCV, 
as  well  as  the  relevance  of  convective  processes  to  the  con¬ 
centration  of  PV  substance,  we  show  PV  horizontal  cross 
sections  at  700  hPa  (Fig.  6),  which  are  taken  from  anima¬ 
tions  that  were  constructed  from  the  MM5  output  on  the  in¬ 
nermost  (highest  resolution)  domain.  The  animation  shows 
the  creation  of  strong,  small  scale,  positive  and  negative  PV 
anomalies  in  the  leading  convective  line  of  the  MCS.  These 
anomalies  are  associated  with  individual  thunderstorm  up¬ 
drafts,  whose  vorticity  has  opposite  signs  due  to  tilting  (Cram 
et  al.,  2002).  After  being  produced  in  the  convective  region, 
the  PV  anomalies  then  move  rearward  into  the  stratiform  pre¬ 
cipitation  region  as  the  convective  line  expands  outward  from 
the  MCV  center  (near  52°  N,  16°  W  in  Fig.  6a).  A  combina¬ 
tion  of  model  diffusion,  vertical  gradients  of  diabatic  heating, 
and  vorticity  center  mergers  results  in  a  smoother  region  of 
positive  PV  emerging  in  the  central  region  of  the  MCV.  For 
the  purpose  of  the  present  discussion,  we  identify  the  central 
region  of  the  MCV  as  the  area  within  the  700  hPa  circulation 
that  is  free  of  strong  negative  PV  anomalies  and  has  a  smooth 
spatial  variation  of  PV  relative  to  the  outer  convective  region. 
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(a) 


(b) 


(c) 


(d) 


Fig.  6.  Trajectory  analysis  and  evolution  of  the  PV  (nr  s-1  Kkg-1)  and  winds  (kt)  on  the  3  km  domain  of  the  MM5  CAPE  simulation  at 
the  following  times  during  the  simulation:  (a)  156. 8h;  (b)  158.4h;  (c)  160.0h;  and  (d)  161. 6h.  The  markings  “If”,  “2f”,  and  “3f”  refer 
to  three  trajectories  originating  within  individual  PV  centers  diagnosed  at  t  =  156.8  h  into  the  simulation,  and  the  curved  arrows  indicate  the 
trajectories  from  f=156.8h  to  the  valid  time  of  the  PV  and  wind  analysis.  Arrows  along  the  trajectories  are  drawn  at  intervals  of  1.6h. 


To  provide  more  quantitative  information  on  the  above 
processes,  we  performed  a  trajectory  analysis  on  three  air 
parcels  originating  from  separate  positive  PV  centers  identi¬ 
fied  in  the  animations.  The  analysis  was  performed  between 
simulation  times  t  =  156.8  h  and  t  =  161.6  h.  The  model  out¬ 
put  interval  for  the  calculation  was  8  min  on  the  innermost 
domain.  To  improve  the  accuracy  of  trajectory  calculations, 
we  used  a  trajectory  time  step  of  80  s.  The  Read,  Interpolate, 
and  Plot  (RIP)  program  (part  of  the  MM5  analysis  system) 
was  used  to  calculate  the  trajectories. 

At  t  =  156.8  h  (Fig.  6a),  the  originating  points  of  these 
three  trajectories,  labeled  “If”,  “2f”,  and  “3f”,  are  shown  (we 
use  this  same  labeling  to  refer  both  to  the  trajectories  and  to 
the  air  parcels  following  them).  By  t  =  158.4  h  (Fig.  6b),  as 
the  convection  expands  outward,  the  three  identified  parcels 
all  initially  move  rearward  relative  to  the  leading  convective 
line  and  begin  to  rotate  cyclonic  ally  about  a  common  center. 


Parcel  3f  ascends  and  moves  east-southeastward,  however, 
eventually  exiting  the  MCS  in  the  strong,  westerly  upper  tro¬ 
pospheric  flow  (Fig.  6d).  At  t  =  160.0h  (Fig.  6c),  parcels  If 
and  2f  continue  to  rotate  about  the  center  of  the  MC  V  and  are 
now  essentially  part  of  the  central  circulation  of  the  MCV.  By 
t  =  161.6  h  (Fig.  6d),  parcel  2f  has  moved  very  close  to  the 
center  of  the  MCV  and  has  begun  to  move  toward  the  south¬ 
west.  Parcel  If  has  rotated  cyclonically  around  the  western 
side  of  the  MCV  center.  The  analyses  presented  in  Figs.  4 
and  5  show  that  the  rotation  comes  about  primarily  due  to  the 
flux  convergence  of  absolute  vorticity  as  the  system  rapidly 
intensifies  due  to  the  diabatic  heating. 

Overall,  the  illustration  shows  that  PV  anomalies  along 
trajectories  If  and  2f  are  created  by  the  leading  edge  con¬ 
vection,  then  move  rearward  and  become  somewhat  less  in¬ 
tense.  It  must  be  noted  that  parcel  2f  descended  from  its 
initial  altitude  near  3000  m  a. g.l.  (Fig.  6a)  to  a  final  altitude 
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Fig.  7.  Back-trajectory  analysis  and  evolution  of  the  PV  (m2  s  1  Kkg  *)  and  winds  (kt)  on  the  3  km  domain  of  the  MM5  CAPE  simulation 
at  the  following  times  during  the  simulation:  (a)  156.8 h;  (b)  158.4h;  (c)  160.0 h;  and  (d)  1 61 .6 h.  The  markings  “lb”,  “2b”,  and  “3b”  refer 
to  three  trajectories  ending  in  individual  PV  centers  diagnosed  at  t  =  161.6  h  into  the  simulation.  The  back-trajectories  were  calculated  from 
f  =  161.6htot  =  156.8  h,  and  the  curved  arrows  indicate  the  same  back-trajectories,  forward,  from  t  =  156.8  h  to  the  valid  time  of  the  PV  and 
wind  analysis.  Arrows  along  the  trajectories  are  drawn  at  intervals  of  1.6 h. 


of  approximately  lOOOma.g.l.,  so  its  final  location  is  be¬ 
low  the  700  hPa  level  shown  in  Fig.  6d.  Parcel  If  remained 
near  3000ma.g.l.,  while  parcel  3f  reached  an  altitude  over 
1 1  000  m  a.g.l.  by  the  end  of  the  analysis  period. 

Since  two  of  the  three  trajectories  ended  at  points  hori¬ 
zontally  or  vertically  displaced  from  the  MCV  central  region 
shown  in  Fig.  6,  we  also  calculated  backward  trajectories 
for  three  points  within  the  700  hPa  MCV  central  region  at 
t  =  161  h  (Fig.  7).  The  three  end  points  chosen  were  from  a 
positive  PV  anomaly  in  the  northern  part  of  the  MCV  (lb), 
a  location  at  the  center  of  the  circulation  (2b),  and  a  positive 
PV  anomaly  from  the  southern  portion  of  the  MCV  (3b).  See 
Fig.  7d  for  the  identified  locations  of  these  three  centers. 

The  backward  trajectory  calculations  show  that  parcel  2b 
and  parcel  3b  were  already  near  the  center  of  the  MCS  at 
the  initial  time  of  t  =  156.8  h  (Fig.  7a).  Both  these  parcels 


rotated  cyclonically  around  each  other  over  the  entire  period 
of  analysis.  Parcel  2b  started  at  an  altitude  of  3600  m  a.g.l., 
and  parcel  3b  started  at  an  altitude  of  3200  m  a.g.l..  Mean¬ 
while,  parcel  lb  started  from  an  altitude  of  2300  m  a.g.l.  at 
a  location  east  of  the  MCS,  entered  the  MCS  in  the  convec¬ 
tive  region  (Fig.  7b),  ascended  to  an  altitude  of  4500  m  a.g.l. 
at  t  =  157.9  h,  then  gradually  descended  to  its  final  altitude 
as  it  became  entrained  into  the  central  region  of  the  MCV 
(Fig.  7d).  All  three  parcels  have  rotated  cyclonically  around 
a  common  center  of  circulation. 

Further  insight  into  the  role  of  stratiform  and  convective 
processes  can  be  gained  by  tracking  the  evolution  of  PV 
along  these  trajectories.  For  this  final  analysis,  we  con¬ 
centrated  only  on  those  trajectories  that  began  (t  =  156.8  h) 
outside  the  MCV  circulation  and  ended  (f  =  161.6  h)  in  the 
central  MCV  circulation  at  an  altitude  between  2500  m  and 
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Time  (h) 

Fig.  8.  Evolution  of  PV  (m2  s  1  Kkg  )  along  parcel  trajectory 
If  (positive  anomaly  identified  in  Fig.  6),  the  trajectory  originating 
from  the  nearest  negative  PV  center  (lfn),  and  the  back-trajectory 
lb  (identified  in  Fig.  7)  according  to  the  simulation  time.  For  lo¬ 
cations  of  these  trajectories  at  various  times,  refer  to  Figs.  6  and  7 
(trajectory  lfn  is  not  shown  but  lies  very  close  to  trajectory  If). 


3500ma.g.l..  Because  our  forward  trajectories  originated 
only  from  large  magnitude  positive  PV  anomalies,  the  anal¬ 
ysis  would  be  biased  to  conclude  that  convection  is  the  pri¬ 
mary  contributor  to  positive  PV  anomalies  in  the  MCV  cen¬ 
tral  region.  Thus,  we  also  calculated,  for  each  of  the  three 
trajectories  in  Fig.  6,  the  corresponding  trajectories  originat¬ 
ing  from  the  nearest  negative  PV  anomaly  (lfn,  2fn,  and  3fn, 
respectively  -  not  shown).  Trajectories  lfn  and  2fn  remained 
very  close  to  trajectories  If  and  2f,  but  only  lfn  finished  be¬ 
tween  2500  m  and  3500ma.g.l..  Thus,  for  our  PV  tracking 
analysis,  we  focused  on  trajectories  lb.  If,  and  lfn. 

The  results  (Fig.  8)  show  that  large  PV  fluctuations  occur 
whenever  the  trajectories  are  located  within  the  MCS.  For 
each  of  the  trajectories,  the  most  rapid  PV  change  occurs  in 
the  earliest  time  intervals,  when  the  parcels  are  located  within 
the  leading  line  of  convection.  Of  particular  interest  is  trajec¬ 
tory  lfn.  PV  values  increase  very  rapidly  through  t  =  157.5  h. 
The  values  then  fluctuate  but  experience  a  net  increase  be¬ 
tween  t  =  157.5  h  and  t  =  161.6  h.  Trajectory  lb  initially  ex¬ 
periences  some  strong  positive  PV  concentration  as  it  enters 
the  convective  region,  but  its  PV  decreases  very  rapidly  over 
a  15  min  time  interval  ending  at  r  =  158  h.  Thereafter,  its  PV 
experiences  a  rather  steady  and  strong  increase. 

The  overall  picture  that  emerges  from  Figs.  5,  6,  7,  and  8 
is  that  the  diabatic  heating  gradients,  and  hence  PV  concen¬ 
tration,  in  the  CAPE  simulation  are  clearly  strongest  in  the 
convective  region  of  the  MCS.  Because  the  tilting  processes 


in  convective  updrafts  generate  dipoles  of  vertical  vorticity, 
the  strong  diabatic  heating  gradients  in  the  convective  region 
are  used  for  the  concentration  of  both  positive  and  negative 
PV  anomalies.  Thus,  it  may  be  difficult  to  conclude,  based 
on  Figs.  6,  7,  and  8,  that  the  ultimate  source  of  positive  PV 
anomalies  in  the  MCV  central  region  is  primarily  from  the 
leading  convective  line.  Regardless,  the  intense  positive  and 
negative  dipoles,  which  require  tilting  processes,  can  be  ex¬ 
pected  to  largely  cancel.  Although  an  exhaustive  attribution 
of  PV  concentration  by  various  processes  is  obviously  a  com¬ 
plex  procedure,  the  largest  source  of  the  positive  PV  anomaly 
in  the  MCV  central  region  has  to  come  from  the  flux  conver¬ 
gence  of  absolute  vorticity  (the  majority  of  which  is  plane¬ 
tary  vorticity  in  the  initial  condition)  acting  in  concert  with 
the  convective  diabatic  heating. 

5  Conclusions 

In  this  paper,  we  have  clarified  several  outstanding  issues 
concerning  the  predominant  mechanism  of  vorticity  gener¬ 
ation  in  mesoscale  convective  vortices  (MCVs)  in  weak  to 
modest  baroclinic  environments  with  nonzero  Coriolis  pa¬ 
rameter.  Using  idealized  mesoscale  numerical  simulations 
of  MCV  development,  we  examined  the  vertical  vorticity 
budgets  in  order  to  quantify  the  contributions  of  flux  con¬ 
vergence  of  absolute  vorticity  versus  tilting  in  the  generation 
of  MCV  vorticity.  We  have  examined  also  the  corresponding 
diabatic  heating  profiles.  By  partitioning  the  diabatic  heating 
between  convective  and  stratiform  regions,  we  elucidated  the 
respective  roles  of  convective  and  stratiform  precipitation  in 
the  concentration  of  potential  vorticity  (PV)  substance. 

The  analyses  indicated  that  the  horizontal  flux  conver¬ 
gence  of  vertical  vorticity  is  the  dominant  mechanism  for 
the  spin-up  and  intensification  of  mid-level  absolute  vortic¬ 
ity.  Correspondingly,  diabatic  heating  and  circulation  bud¬ 
gets  demonstrated  that  the  vertical  gradient  of  diabatic  heat¬ 
ing  is  supportive  of  low-  to  mid-level  concentration  of  PV 
substance.  When  the  CAPE  is  representative  of  warm  season 
conditions,  and  during  the  early  stages  of  MCV  development, 
convective  precipitation  plays  the  dominant  role  in  the  con¬ 
centration  of  PV.  During  later  stages  in  the  MCV  life  cycle, 
the  contribution  from  stratiform  processes  increases.  In  con¬ 
trast,  in  an  environment  with  near  zero  CAPE,  stratiform  pre¬ 
cipitation  is  the  dominant  contributor  to  PV  concentration. 

Trajectory  analyses  and  plots  of  PV  on  the  innermost  grid 
of  the  simulations  revealed  that  the  convective  region  has 
strong  PV  anomalies  of  positive  and  negative  sign.  If  we 
use  the  absolute  vorticity  diagnostics  as  a  proxy  for  the  evo¬ 
lution  of  PV  substance,  then  we  find  a  cyclonic  vorticity  rich 
environment  favors  the  emergence  of  a  cyclonic  PV  anomaly 
behind  the  leading  convective  line. 
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